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A possible experimental observation of the chiral magnetic effect in heavy ion collisions at RHIC 
was recently reported by the STAR Collaboration. We study signatures of this effect in SU(2) 
lattice gluodynamics with the chirally invariant Dirac operator. We find that at zero temperature 
\Q the local fluctuations of an electric current of quarks and chirality fluctuations increase with exter- 



nal Abelian magnetic field. The external magnetic field leads to spatial separation of the quark's 
electric charges. The separation increases with the strength of the magnetic field. As tempera- 



Q\ ture gets higher the dependence of these quantities on the strength of the magnetic field becomes 

weaker. In the deconfinement phase the local fluctuations of the chiral density and of the spatial 
components of the quarks electric current are large and are almost independent on the external 
magnetic field. The local fluctuations of the electric charge density decrease with the strength of 
£3 the magnetic field in this phase. 
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1. Introduction 

Very strong magnetic fields of the hadronic scale can significantly modify the properties of 
strongly interacting matter: they change the order of the phase transition from the confinement 
phase to the quark-gluon plasma, shift the position of the transition line [1], etc. At the Relativistic 
Heavy Ion Collider (RHIC) in noncentral heavy-ion collisions the strong magnetic field arises due 
to the relative motion of the ions and products of the collision [2]. The induced magnetic field is 
perpendicular to the reaction plane. At first moments (t ~ 1 fm/c) of the collision the value of 
magnetic field at RHIC may reach the hadronic scale, ^fqB ~ (10 — 100 MeV) [2, 3]. Such strong 
magnetic fields can also be created in future at the ALICE experiment at LHC, at the Facility for 
Antiproton and Ion Research (FAIR) at GSI, and at the Nuclotron Ion Collider fAcility (NICA) in 
Dubna. 

The strong magnetic fields can also enhance the chiral symmetry breaking. Chiral perturbation 
theory predicts a linear rise of the chiral condensate (yy) with the strength of magnetic field [4], 
that is in accordance with our lattice results [5]. In the AdS/QCD approach [6], in Nambu-Jona- 
Lasinio model [7] and in sigma models [8] the value of (Y¥) rises quadratically with the external 
magnetic field. 

Strong magnetic fields lead also to the chiral magnetization of QCD vacuum at zero and finite 
temperature. This effect has a paramagnetic nature because the quarks are 1/2-spin particles and 
the external magnetic field leads to polarization of magnetic moments of quarks in the direction 
of the magnetic field. The magnetization of the QCD vacuum is essential for the properties of 
the nucleon magnetic moments [9] and other nonperturbative features of hadrons [10]. Our lattice 
simulations [11] demonstrate at the relatively weak magnetic fields the chiral magnetization rises 
linearly in accordance with Ref. [9]. At larger values of the magnetic fields the magnetization 
depends nonlinearly on the magnetic field [12]. 

We also observed an anomalous quark magnetization in intense magnetic fields which appears 
due to the topological fluctuations of the QCD vacuum [13]. The essence of this effect is that 
quark gets a local electric dipole moment in addition to the magnetic moment. Both magnetic and 
anomalous electric moments are parallel to the axis of the external magnetic field. The generation 
of the anomalous electric dipole moment is a spin analogue of the chiral magnetic effect (CME). 

CME is the generation of a local electric current in the direction of the external magnetic field 
in topologically nontrivial configurations of the gauge fields [2, 14]. If we consider u and J-quarks 
as massless particles then the right-handed quarks should move in the direction of the magnetic 
field and the left-handed quarks should move in the opposite direction because in the external field 
magnetic moments of quarks are parallel to the direction of the field. Nonzero topological charge 
of gauge fields leads to a local imbalance between left-handed and right-handed quarks, which, in 
turn, gives rise to a nonzero net electric current along the axis of the magnetic field. An evidence 
of the CME was observed by the STAR Collaboration at RHIC as an non-statistical asymmetry 
of negatively and positively charged particles emitted at different sides of the reaction plane of 
noncentral heavy-ion collisions [15, 16]. The observed effect indicates the presence of the chirality 
fluctuations in heavy-ion collisions. Below we report the evidence of the CME in the lattice gauge 
theory following our original study [17]. 
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Figure 1: The excess of the electric charge density due to the external magnetic field qB — 0.7 GeV 2 (left) 
and qB = 1.81 GeV 2 (right) on 14 4 lattice corresponding to the physical volume of (1.44 fm) 4 . The red 
regions corresponds to the excess of the positive charge density, while the blue color marks the excess of the 
negative charge density. 

2. Details of calculations 

We generate statistically independent SU(2) gauge field configurations with the tadpole im- 
proved Symanzik action [18]. Using the chirally invariant Dirac operator [19], we solve the Dirac 
equation Dy/& = X^Yk numerically and determine the corresponding eigenf unctions \j/k and eigen- 
values X^. Here D = y^(d^— iA^ ) is the massless Dirac operator in the gauge field . The uniform 
magnetic field is introduced into this operator as described in Ref. [5]. 

We performed the zero-temperature simulations on 14 4 lattice with the lattice spacing a = 
0.103 fm and on 16 4 lattice with the lattice spacings and a = 0.103 fm and a = 0.089 fm. At finite 
temperature we used 16 3 x 6 lattices with the spacings a = 0.128 fm (T = 256 MeV = 0.82 T c ) 
and a = 0.095 fm (T = 350 MeV = 1.12 T c ). The critical temperature of the phase transition in 
SU (2) gauge theory is T c = 313.(3) MeV [20]. The values of the magnetic field are quantized due 
to the periodic boundary conditions imposed in finite lattice volume. In our simulations a minimal 
nonzero value of magnetic field is qB m i n = (348 MeV) 2 . 

3. Fluctuations of electric charge density 

In order to study the simplest signatures of the chiral magnetic effect we first explore the 
fluctuations of the electromagnetic current 



The zeroth component of (3.1) corresponds to a local charge density. We remove the diverging 
ultraviolet contributions by subtracting the charge density at zero magnetic field B from the charge 
density calculated at nonzero B. The result of the subtraction is a nonperturbative infrared quantity 
which is plotted in Fig. 1. Here we visualize the level surfaces of two typical charge density 
distributions in a fixed time slice of the zero-temperature gauge field configuration for two values 
of the external magnetic field B. The magnetic field is directed from the bottom to the top. 



JmW = W)7iiW{x). 



(3.1) 
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Figure 2: The expectation values of the squared chirality (left) and of the squares of the longitudinal (jo and 
73 ) and the transverse (J\ and j'2) electric currents (right) vs the external magnetic field at three temperatures. 

We see in Fig. 1 that the local charge density grows with the external magnetic field. The 
regions with larger charge density are extended along the direction of the magnetic field because at 
a sufficiently strong field the quarks occupy the lowest Landau level. In the absence of the gluon 
background the trajectories of the quarks should be elongated along the direction of the magnetic 
field. The gluon field distorts the motion of quarks, so that the trajectories of quarks diffuse. 

4. Fluctuations of chirality 

The other quantity which characterizes the chiral magnetic effect is the local chirality: 

p 5 (x) = \jf(x)Y5V(x). (4.1) 
The average value of (4.1) equals to zero, therefore we measure the squared chirality: 

(pI)ir(B,T) = ^ / d 4 x(p 5 (x) 2 ) BJ - i / d 4 x( P5 (x) 2 ) B=0iT = , (4.2) 
V Jv V Jv 

where V is the lattice volume. The subtraction of the B = and T = quantity is performed in 
order to regularize the total result in the ultraviolet region. 

Fig. 2 (left) shows the squared chirality at different temperatures T. At zero temperature the 
chirality fluctuations increase quickly with the value of B. In the confinement phase at nonzero 
temperature (T = 0.82 T c ) the grows rate becomes smaller and in the deconfinement phase (T = 
1.12 T c ) the squared chirality is almost independent on the value of the magnetic field. 

The value of the chirality fluctuations at zero temperature and sufficiently large magnetic fields 
(qB ~ lGeV 2 ) is close to the value of (p|) in the deconfinement phase (T = 1.12 T c ) at zero mag- 
netic field. Thus, we can conclude that strong chirality fluctuations and hence the chiral magnetic 
effect could be also observed in a cold nuclear matter exposed to the strong magnetic fields. How- 
ever, in this case the experimental signatures of the CME can be quite different from those in 
heavy-ion collisions. 
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Figure 3: The squares of the longitudinal and transverse components of the electric currents vs the external 
magnetic field in the confinement phase (left) and in the deconfinement phase (right). 



5. Fluctuations of electromagnetic current 

We study the fluctuations of the electric current: 

(jl)m(B,T) = i J d 4 x(jfa)) B ,T -^J d 4 x(jl(x)) B=0 ,T=o- (5.1) 

We checked on the lattice, that all components of the electric current (3.1) are on average equal to 
zero within statistical error bars. The component of the current along the direction of the magnetic 
field is called longitudinal. Transverse components are perpendicular to the direction of the field. 
At T = the magnetic field in the fL = 3 direction breaks the rotational symmetry, so that we have: 
(jj) = (jj) and (f ) = (jj). 

In Fig. 2 (right) we show all components of the electric field. All of them grow with the 
strength of the field. The longitudinal component of the electric current and the electric charge 
fluctuate stronger than the transverse ones. The fluctuations of the spatially transverse components 
grow with the field because the transverse momentum of a quark occupying the lowest Landau 
level increases with the strength of the magnetic field. 

At T the rotational symmetry is broken by the temperature and by the external magnetic 
field, so that the longitudinal and zero components are no more equivalent, (j'q) ^ while the 
spatially transverse components are still degenerate, (jf) = (fy. 

In Fig. 3 (left) we show the fluctuations of the electric current in the confinement phase. At 
weak fields all components decrease slightly (the decrease is, however, within the statistical errors). 
Then the fluctuations start to increase with the value of the magnetic field. The components (j'q) 
and (jg) are larger than and (jf). From Fig. 3 (right) we see, that in the deconfinement phase 
the components (j^) and (jj) are almost independent on the external magnetic field B, while 
the charge fluctuation (Jq) is a decreasing function of B. This observation can be explained by the 
Debye screening in the deconfinement phase. 
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Figure 4: Left: the correlation function of p| (x) and j^(x), Eq. (6.1). Right: The values of (p|) and (jfy at 
an instanton-like configuration with 2=1. 

6. Chirality-current correlations 

We observed that the (psjp) = 0. This fact is natural due to the CP-oddness of this quantity. 
Next, we study the correlation between the squares of the chirality and the induced current: 



_ (pljl)-(pl )Ul 



<lPFJi) • • (6-1) 



We found that in the confinement phase the value of (6.1) vanishes for all components of the 
electric currents. In the deconfinement phase the correlator is decreasing function of the field 
[Fig. 4 (left)]. Therefore, the enhancement of the fluctuations of the electric current in the direction 
of the field is not locally correlated with the chirality fluctuations. 

7. Instanton-like configuration 

We have also studied signatures of the chiral magnetic effect at an instanton-like gauge field 
configuration with unit topological charge, 2=1. In Fig. 4 (right) we show the fluctuations of 
the local chirality and the fluctuations of the squared longitudinal electric current in the instanton 
background. We observed the growth of the chiral density with the strength of the applied external 
magnetic field. The grows is caused by a difference between the particles with positive and negative 
chiralities in the Q = 1 background. The net longitudinal electric current, jl = j'q + j|, increases 
with the strength of the field as well. This is the signature of the chiral magnetic effect. 

8. Conclusion 

We found the signatures of the chiral magnetic effect in lattice SU(2) gluodynamics in a back- 
ground of the strong external magnetic field. This effect arises due to imbalance between left- and 
right-handed quarks (the local chirality). The average chirality in a large enough volume is equal 
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to zero, but the local fluctuations of the chiral density can be sufficiently strong. We found that the 
fluctuations of the chirality at zero temperature increase with the strength of the magnetic field. As 
the temperature increases the fluctuations become weaker dependent on the strength of the field B. 
In the deconfinement phase the local chirality fluctuation are almost independent on B. 

We found that in the confinement phase all components of the electric current of quarks are 
growing in intense magnetic fields. The longitudinal (with respect to the axis of the magnetic 
field) currents fluctuate stronger than the transverse ones. In the deconfinement phase the thermal 
fluctuations lead to the increase of the current fluctuations while all spatial components of the 
current are insensitive to the value of the magnetic field. The charge fluctuations decrease as the 
field gets stronger. The fluctuations of the electric current and chirality are almost independent. 

Finally, we observed the chiral magnetic effect at the instanton-like configuration. 
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